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ABSTRACT 


The mini-proceedings of the 17*^ Meeting of the ’’Working Group on Radiative Corrections 
and MonteCarlo Generators for Low Energies” held in Frascati, 20**^ - 21®* April, are pre¬ 
sented. These meetings, started in 2006, have as aim to bring together experimentalists and 
theoreticians working in the fields of meson transition form factors, hadronic contributions 
to the anomalous magnetic moment of the leptons, and the effective fine structure constant. 
The development of MonteCarlo generators and Radiative Corrections for precision 
and r-lepton physics are also covered, with emphasis on meson production. Heavy quark 
masses were covered as well in this edition. 

The web page of the conference: 

https://agenda.infn.it/conferenceDisplay.py?ovw=True&confId=9287 
contains the presentations. 

We acknowledge the support and hospitality of the Laboratori Nazionali di Frascati. 
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1 Introduction to the 17^^ Radio MontecarLow Work¬ 
ing Group meeting 


H. Czyz^ and G. Venanzoni^ 

^Institute of Physics, University of Silesia, 40007 Katowice, Poland 
^Laboratori Nazionali di Frascati dellINFN, 00044 Frascati, Italy 

The importance of continuous and close collaboration between the experimental and 
theoretical groups is crucial in the quest for precision in hadronic physics. This is the 
reason why the Working Group on “Radiative Corrections and Monte Carlo Generators for 
Low Energies” (Radio MonteCarLow) was formed a few years ago bringing together experts 
(theorists and experimentalists) working in the held of low-energy e’''e“ physics and partly 
also the r community. Its main motivation was to understand the status and the precision of 
the Monte Carlo generators (MC) used to analyze the hadronic cross section measurements 
obtained as well with energy scans as with radiative return, to determine luminosities, and 
whatever possible to perform tuned comparisons, i.e. comparisons of MC generators with a 
common set of input parameters and experimental cuts. This main effort was summarized 
in a report published in 2010 [T]. During the years the WG structure has been enriched of 
more physics items and now it includes seven subgroups: Luminosity, R-measurement, ISR, 
Hadronic VP g — 2 and Delta alpha, gamma-gamma physics, FSR models, tan decays. 

During the workshop the last achievements of each subgroups have been presented. The 
present accuracy and the future prospects of MC generators for e+e“ into leptonic, 77 , 
hadronic and tan hnal states have been reviewed. Recent proposals on the hadronic LO 
contribution to the g — 2 oi the muon in the space-like region and on a dispersive formalism 
for 7*7 —)■ 7171 have been discussed. A new method for extracting heavy quark masses 
from BES and CELLO data was as well discussed. New results from CMD3 and BESIII 
experiments have been presented, together with new Monte Carlo generator developments. 

The present workshop, being its 17*^ edition, was held from the 20*^ to the 21 ^* April 
2015, at the Laboratori Nazionali di Frascati dellINFN, Italy. 

Webpage of the conference is 

https://agenda.infn.it/conferenceDisplay.py?ovw=True&confId=9287 

where detailed program and talks can be found. 

All the information on the WG can be found at the web page: 

http://www.lnf.infn.it/wg/sighad/ 
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Summaries of the talks 


2.1 Recent Results from VEPP-2000: Data and Generators 

S. Eidelman 

Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090, Russia, 
Novosibirsk State University, Novosibirsk 630090, Russia 


Two detectors at the VEPP-2000 e+e collider in Novosibirsk, CMD-3 and SND, con¬ 
tinued data processing of an integrated luminosity of 60 fb ^ collected by each in 2011 — 
2013 in the center-of-mass (c.m.) energy range 320 — 2000 MeV. Lately analysis has been 
mainly focused on the c.m. energies above the 0 meson. In particular, both groups made 
an attempt to improve the precision of various cross sections important for the muon g — 2 
problem pQ. 

Measurements of an integrated luminosity at CMD-3 use two processes, e+e“ e~^e~ 
and e’''e“ ^ 77 , allowing a precision of ~ 1% [2]. At SND, events of large-angle Bhabha 
scattering are used to determine an integrated luminosity with a systematic accuracy of 

2% |3]. 

Both detectors continued analysis of the collected data samples and reported measure¬ 
ments of cross sections for various processes with pions and 7 mesons. SND published hnal 
results on e^e“ —?■ 7 r+ 7 r“ 7 r‘’ [1] and tt+tt”?] [3]. CMD-3 is close to obtaining hnal 

values of cr(27r+27r“27r°) [5]. Analysis of dynamics should be performed simultaneously with 
a study of the 37 r+ 37 r“ hnal state, which cross section had been reported before [ 6 ]. 

CMD-3 continued studies of various processes with kaons in the hnal state using good 
K/n separation based on measuring dE/dx in the drift chamber. In addition to the process 
e’''e“ K^K~ [7], they reported close to hnal results on the cross section and dynamics 
of the and K^K~r] hnal states. They have also presented their hrst results on 

the process e’''e“ —)■ 

SND has already published their measurements of the processes with only neutral par¬ 
ticles in the hnal state - —)■ [S] and e+e“ —)■ 77 [H]- 

A study of the nucleon form factors near threshold was continued. SND signihcantly im¬ 
proved the precision of a{e'^e~ —?■ nn) [10] compared to the previous results from FENICE m 
CMD-3 measured the cross section of the process e+e —)■ pp and made an attempt to ex¬ 
tract the ratio of the electric and magnetic form factors based on the angular distribution 
of the hnal nucleons [T2] . 

Both detectors used an original method of Ref. [TB] to measure the partial width of a 
strongly suppressed p' —)■ e+e“ decay using the inverse process. CMD-3 reported an upper 
limit of T{p' —)■ < 0.0024 eV at 90% CL based on 2.69 pb“^ and one mode of p' 

decay [H]. SND used 2.9 pb“^ and hve modes of p' decay to improve it to < 0.0020 eV. 
Finally, they combine the data samples of CMD-3 and SND to hnd r(? 7 ' —)■ < 0.0011 

eV at 90% CL [T3] which is still about two orders of magnitude below the unitary bound. 
SND has also performed a feasibility study for a search for p e^e~ via e^e“ p and 
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concluded that the only promising decay mode for that is —)■ A dedicated two-week 

run with the luminosity expected at the c.m.energy around the rj meson mass will allow to 
improve the existing limit [TB] . 

A need of high-precision measurements of hadronic cross sections and studies of rich 
resonance dynamics demand Monte Carlo (MC) generators, which properly take into account 
interference effects and all symmetries. Progress in development of the generic MC generator 
of —)■ hadrons used by CMD-3 for background studies was described na. A generator 

for production of three pseudoscalars in e’''e“ annihilation found further applications [TB] . 
Intensively discussed was an important question about radiative corrections included in 
these and similar MC generators. 
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2.2 Improving the luminosity measurement at BESIII 
using the Bhabha event generator Babayaga@NLO 

A. Hafner 

Institut fiir Kernphysik, Johannes Gutenberg Universitat Mainz, 
J.-J.-Becher-Weg 45, 55099 Mainz, Deutschland. 


There is a long standing discrepancy between the experimental measurement of the 
anomalous magnetic moment of the muon PP and its theoretical prediction within 
the Standard Model of particle physics. This discrepancy amounts to 3-4 standard devia¬ 
tions [SIE]. New experimental measurements BE] will reduce the experimental uncertainty 
by approximately a factor of 4 within this decade. The dominant uncertainty on the the¬ 
oretical prediction of the anomalous magnetic moment of the muon stems from the 
contribution of hadronic Vacuum Polarization (VP), = (692.3 ± 4.2) • 10“^° [2]. From 

causality and analyticity of the VP amplitude a dispersion relation for the VP contribution 
to can be derived [6]. This relation requires the inclusive hadronic cross section as 
input. The largest weight is given to low energy contributions. Therefore, the reaction 
a{e'^e~ —)■ tt+tt”) contributes with approximately 75%. 

The standard experimental approach is to measure the required hadronic cross sections 
exclusively at e^e“ energy scan experiments. Since the last decade, the method of Ini¬ 
tial State Radiation (ISR) is used as an alternative approach to measure cross sections of 
exclusive final states at high luminosity flavor factories, running at a hxed center-of-mass 
energy. The emitance of a high energy photon from initial state opens the window to low 
energy hadron physics. KLOE, running on the 0 resonance, measured the —)■ tt+tt” 

hnal state [7] with a precision of better than 1% in the peak region. BABAR, running on 
the T(4S') resonance, has an extensive ISR-scan program with various hnal states up to six 
hadrons from energy threshold up to 4.5 GeV [8]. The BABAR measurement with a precision 
of 0.5% of the TT+TT” hnal state shows a discrepancy in and above the p region of up to 
2-3 standard deviations to the KLOE measurement. Due to this diherence, the resulting 
uncertainty for is similar to the uncertainties of the individual measurements. 

Energy scan measurements of the 7r'''7r“ cross section by GMD-3 and SND experiments 
are expected in the near future with an aimed uncertainty of below 1% and 0.5% respec¬ 
tively. In addition, a new ISR measurement at BESIII, running in the charmonium region 
and below, is currently performed. The data is statistically and systematically competitive 
to BABAR and KLOE. In a hrst step, the aim is a precision of 1% in the p-resonance region 
using a data sample at the -^(3770). 

The dominating source of systematic uncertainty of this cross section at BESIII stems 
from the uncertainty of the luminosity measurement, which at present by itself amounts to 
1 % M- The precision was achieved by measuring Bhabha events and using theoretical input 
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of event generators. Thus, in this measurement, the main contributions to the uncertainty 
stem from the uncertainty due to the event generator Babayaga. 3.6 Hi of 0.5% and a very 
conservative approach for the uncertainty due to the polar angular requirement of 0.7%. 
The uncertainty of the event generator can be reduced to 0.1% by using the successor of 
the event generator, Babayaga@NLO m- First studies of the polar angular distribution 
clearly indicate, that the uncertainty due to the corresponding requirement can also be sig- 
nihcantly reduced. Effects due to VP uncertainties on the Bhabha cross section, and thus 
the luminosity measurement, have been taken into account and are shown to be negligible 
for the rather broad '0(3770) resonance. 

With Babayaga@NLO and the additional systematic studies, the uncertainty for the 
luminosity will be signihcantly reduced, allowing the cross section measurement cr(e’''e“ 
7r~^7r~) at BESIII to reach a precision of 1% and below. 

I want to thank the organization committee of the 17*^ Radio Monte Carlo working group 
meeting for the warm hospitality in Frascati and the participants for the fruitful discussions. 
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2.3 Measurement of the e+e —tt+tt cross section 
using initial state radiation at BESIII 

C. F. Redmer , B. Kloss, and A. Denig 

Institut fiir Kernphysik, Johannes Gutenberg Universitat Mainz, 
J.-J.-Becher-Weg 45, 55099 Mainz, Deutschland. 


Precise measurements of hadronic cross sections at e’''e“ colliders are an important input 
to the Standard Model prediction of the muon anomaly a^. Currently, the accuracy of this 
prediction is entirely limited by the understanding of the hadronic contributions to pQ. 
The largest of these contributions stems from the hadronic vacuum polarization contribu¬ 
tion. It is handled within a dispersive framework, which requires a{e'^e~ hadrons) as 
experimental input [2]. The cross section airir = cr{e'^e~ —)■ 7r’''7r“) contributes more than 
70% to this dispersion relation and is, hence, by far the most important exclusive hadronic 
channel in the endeavor of improving the knowledge on the hadronic contributions to a^. 

To date, the two most accurate measurements of have been obtained by the KLOE [3], 
and the BABAR |1] collaborations. Both experiments claim a precision of better than 1% 
in the energy range below 1 GeV, however, a discrepancy of approximately 3% on the peak 
region of the p{770) resonance is observed. The discrepancy is even increasing towards 
higher energies and has a large impact on the SM prediction of [2]. Thus, an indepen¬ 
dent measurement of competitive accuracy, i.e. in the order of 1 %, is required to settle the 
issue. 

This measurement has been performed at the BESIII experiment [5], operated at the 
symmetric e’''e“ collider BEPGII in Beijing, Ghina. Based on 2.9 fb“^ of data [ 6 ], taken 
at a center-of-mass energy y/s = 3.773 GeV, the method of initial state radiation (ISR) 
has been exploited, measuring events of the type Tr+vr'y. In this way, the mass 

range between 600 and 900 MeV/c^, which corresponds to the important p peak region, was 
studied. 

The production of muon pairs constitutes the main background contribution in this 
analysis. To suppress events of the kind an artihcial neural network [7] has 

been trained and tested. Monte Garlo (MG) simulations of the signal and background final 
states, 7 r’'' 7 r “7 and /i+/i“ 7 , have been used as input. The Phokhara generator [ 8 ] has been 
used to generate signal and background samples. Possible discrepancies between data and 
MG simulation due to imperfections in the detector simulation, have to be considered. To 
accomplish a high precision analysis, track-based data-MG correction factors were obtained, 
comparing the analysis of nearly background free pion and muon samples in data with MG 
simulations. 

Two independent normalization methods can be used to extract On the one hand, 
the efficiency corrected number of Tr+vr^y events can be normalized to the luminosity 
and the radiator function [9]. On the other hand, the R ratio can be determined, i.e. is 
normalized to the number of p~^p~'y events. Both methods have been used in the analysis, 
and agree within the errors. However, the the final result has been obtain based on the first 
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method, since the second one is limited by the statistics. 

The extracted cross section can be used to calculate the two-pion contribution to a^. 
The preliminary result in the mass range between 600 and 900 MeV is found to be = 

374.4 ± 2.6stat ± 4.9<jj;s- Compared to the values published by the BaBar and KLOE collab¬ 
orations for the same energy range, the BESIII result almost coincides with their average 
value and agrees with both results within errors. The preliminary estimated systematic 
uncertainty of 1.3% of the BESIII measurement is dominated by the uncertainty of the 
luminosity measurement of 1% [B]. We are conhdent to reduce this uncertainty by thor¬ 
oughly re-evaluating the luminosity measurement [10]. Nevertheless, this result conhrms the 
discrepancy in between SM and experiment on the level of 3 to 4 standard deviations. 
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2.4 Primary Monte-Carlo generator of the process e'^e — /o(1370)/9(770) 
for the CMD-3 experiment 

P.A. Lukin^, V.E. Lyubovitskij^ 

^ Budker Institute of Nuclear Physics and Novosibirsk State University 
^ Tuebingen University, Tomsk State University, Tomsk Polytechnic University 


Electron-positron collider VEPP-2000 [T] has been operating in Budker Institute of Nu¬ 
clear Physics since 2010. Center-of-mass energy range covered by the collider is 

from threshold of hadron production and up to 2 GeV. Special optics, so called “round 
beams”, used in the collider construction, allowed to obtain luminosity 2 x 10^^ cm“^-s“^ at 
Ec.m. =1.8 GeV. 

The general purpose detector GMD-3 has been described in detail elsewhere [2]. Its 
tracking system consists of a cylindrical drift chamber (DG) |3] and double-layer multiwire 
proportional Z-chamber, both also used for a trigger, and both inside a thin (0.2 Xq) su¬ 
perconducting solenoid with a held of 1.3 T. The liquid xenon (LXe) barrel calorimeter 
with 5.4 Xq thickness has hne electrode structure, providing 1-2 mm spatial resolution [1], 
and shares the cryostat vacuum volume with the superconducting solenoid. The barrel Gsl 
crystal calorimeter with thickness of 8.1 Xq is placed outside the LXe calorimeter, and the 
end-cap BGO calorimeter with a thickness of 13.4 Xq is placed inside the solenoid [Sj. The 
luminosity is measured using events of Bhabha scattering at large angles [B] . 

Physics program of the GMD-3 experiment includes the study of the multi-hadron 
production. The cross section measurement of the e+e“ 3(7r+7r“) process in Ec.m. = 

1.5-2.0 GeV has been already published [7]. Preliminary results the study 2(7r’''7r“7r°) 

hnal state has been reported [S]. 

The study of intermediate states which lead to 2{7r~^7r~7r^) hnal state is essential to 
correctly describe the angular correlations between the particles and determine the regis¬ 
tration efficiency of the process under study. As it was reported at [9] the intermediate 
states Ci;(782)37r, ci;(782)?7(545) and p(770)(47r)5_^a.ue allow satisfactorily describe mass and 
angular distributions of the 2(7r+7r“7r°) production in Ec.m. = 1-5-1.7 GeV. 

For higher Ec.m. h was found possible to describe 7^(545) signal, seen in three-pion mass 
distribution of the 2(7r'''7r“7r°), by the process —)■ ao(980)p(770) with dominant decay 

of ao(980) into 17 ( 545)77 as it is described in [10]. The primary Monte-Garlo generator for 
the process has been created out and implemented into the GMD-3 experiment Monte-Garlo 
simulation package. Using the generator the signal of the e+e“ —)■ ao(980)p(770) process 
has been observed in the experimental data for 2{7r~^7r~7r^) hnal state. 

In the present study, it was made an attempt to describe 27r- dvr and dvr mass distributions 
as well as angular correlations for 2(7r’''7r“7r‘^) hnal state at Ec.m. > 1-8 GeV by contributions 
of Ci;(782)37r, ao(980)p(770) and /o(1370)p(770) intermediate states. Two decay channels 
of /o(1370) have been studied /o —)■ p’''(770)p“(770) and /o(1370) —?■ 7r+7r“27r°. It was 
found that process —?■ /o(1370)p°(770) with subsequent decay /o p'^{770)p~{770) 

reasonably describe both masses and angular distributions of the —)■ 2(7r+7r“7r°) process 
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at Ec.m. = 2.0 GeV. And we could not described masses and angular distributions of the 
—)■ 2(7r+7r“7r°) process at A'c.m. = 1-8 GeV by neither /o(1370) p+(770)p“(770) nor 

by /o(1370) 7r+7r“27r° decays. 
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2.5 Event Generators at Belle II 


T. Ferber 

DESY, D-22607 Hamburg, Germany 


The next generation B-factory Belle II at the upgraded KEKB accelerator, SuperKEKB, 
is aiming to start data taking in 2017. The broad physics program covers e.g. physics with 
B and D mesons, /i and r leptons as well as measurements using the method of radiative 
returns and direct searches for new physics. The expected dataset will exceed the one col¬ 
lected by the predecessor Belle by a factor of 50 and imposes high precision requirements 
on the used event generators and requires a flexible and powerful software framework. The 
bulk of the data will be collected at the T (dS*) resonance, but it is planned to collect sizable 
data sets also off-resonance and at energies around the narrow resonances T(IS'), T(2S') 
and T(3^). 

Belle II uses a single software framework, basf2 HlEl, for all data processing tasks which 
runs on standard Linux systems. It is based on a user-dehned chain of individual modules 
where each subsequent module can read data from the preceding modules from a so-called 
data store. Event generators usually serve as specialized modules that provide four vectors 
to be fed into the subsequent GEANT4 based detector simulation as well as precision cross 
section calculations used for normalization. If available from the generator, all mother- 
daughter relations of both unstable particles and radiative photons are stored and used in 
Monte Garlo truth matching during analysis. FORTRAN-based generators are interfaced 
using extern “G” functions where the user inputs and generator outputs are provided as 
FORTRAN common blocks and global G/G-I--I- extern structs of the same name. All gener¬ 
ators use a random generator provided by the basf2 framework. The original input interfaces 
are replaced by Python steering scripts that provide access to the generator options in a 
uniform and user-friendly way. 

The available physics generators include BABAYAGA.NLO E , BHWIDE |1], KKMG4.19 [5l 
E], PHOKHARAQ.lb [7], KORALW1.51 [8], AAFH |9], BBRE^ and EvtGen [H]. 

TAUOLA [12] and PHOTOS [13] are used by EvtGen and KKMG to handle r decays and 
radiative corrections in decays. MadEvent [11] is used to simulate New Physics processes. 
Light quark continuum is modeled using KKMG (hard interaction), PYTHIAS (fragmenta¬ 
tion) [15] and EvtGen (decays). Standardized HepEvt or Les Houches event (LHE) format 
can be read by dedicated input modules. The data from the data store can be obtained at 
any stage of the module-chain in ROOT or HepEvt hie format. 

In conclusion, all basic event generators are available in basf2 and ready for physics and 
trigger studies. Future projects will focus on the precision validation and improvements 
needed to match the demanding precision requirements for luminosity measurements and 

^C-|--|- implementation based on the original FORTRAN code m- 


13 




low multiplicity physics especially at the narrow resonances. A semi-automatic framework 
to check EvtGen models is under development. Work has started to use the Belle datasets 
taken off-resonance and at the T(IS') to tune PYTHIAS within the basf2 framework. 
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2.6 Heavy quark masses form QCD Sum Rules 


P. Masjuan 


PRISMA Cluster of Excellence, Institut fiir Kernpliysik and Helmholtz Institut Mainz, 
Johannes Gutenberg-Universitat Mainz, D-55099 Mainz, Germany 

Heavy quark masses are interesting not only for being fundamental parameters not pre¬ 
dictable within the Standard Model (SM), but for their implications on many phenomeno¬ 
logical scenarios as well. Being a free parameter in the SM, their values should be extracted 
from experimental information which already demands an important synergy between both 
experiment and theory. Beyond that, the heavy quark mass enters into decays, from Higgs 
(where the mass enters squared) to B decays (where enters with the 5th exponential). For 
precision physics, an accurate knowledge of that important input is mandatory, specially if 
one foresees Higgs branching ratios at the per mil level. 

From the different methods to extract the quark mass from experimental data, we follow 
in Ref. pj the Sum Rule s one [a [31 HE]. The method relates, thanks to the optical theorem, 
the transverse part of the correlator of two heavy-quark vector currents denoted by nq(t) 
—where the caret indicates MS subtraction— with integrals over the measured ratio R, the 
ratio between the —?■ hadrons) over the a{e^e~ —)■ 

The nij(t) can be calculated order by order within perturbative QCD as an expansion 
of the strong coupling a(s), and obeys the subtracted dispersion relation [6] 



( 1 ) 


where Rq{s) = 127rlmflq(s), and is the mass of the heavy quark that we want to deter¬ 
mine. In the limit t —)■ 0, Eq. ([T]) coincides with the hrst moment, AJi, of Ilg(t). In general, 
there is a sum rule for each higher moment as well laEiiiiisi: 



Taking the opposite limit, i.e. t oo, of Eq. ([T]) after multiplying with t and with a proper 
regularization, one can dehne a sum rule for the 0*^ moment [Bj. At a given, hxed order 
of pQCD, the required regularization can be obtained by subtracting the zero-mass limit of 
Rq, which we denote 3QqAi(s) with Qq the quark charge. A^(s) is known up to 0{al) [7|. 

By the optical theorem, Rq can be related to the measurable cross section for heavy- 
quark production in annihilation. We assume that below the threshold, the cross 

section is determined by a small number of narrow Breit-Wigner resonances [2] 



( 3 ) 
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The masses Mr and electronic widths T^j of the resonances are collected from the PDG [8] 
and aem(^ij) is the running hne structure constant at the resonance. Then, Rq{s) = 
with accounting for the continuum production region. 

Invoking global quark-hadron duality, we also assume that continuum production can 
be described on average by the simple ansatz [6] (which coincides with the reconstruction 
of the correlator at the given order piiniiii]) 


= 3QlX\{s) 



Am^{2M) 

s' 


1 + 


(4) 


where s' = s -|- 4(mg(2M) — M^) and Trig is the running MS heavy-quark mass, evaluated 
at 0{ag) [12] at the scale 2M, M taken as mass of the lightest pseudoscalar heavy meson. 
Ag is a constant to be determined. interpolates smoothly between the threshold 

and the onset of open heavy-quark pair production. It coincides asymptotically with the 
prediction of pQCD for massive quarks. Two parameters are free to be determine after hts 
to the experimental data by BES [131 HU [131 HSl HZ], and CLEO [13] on the open charm 
region: Ag and ihg, for which two different sum rules will be used, the 0th and the 2nd. In 
Ref. [I], not only the central values of the heavy masses are given, but also an exhaustive 
scrutiny of the many error sources is also presented. 

In this talk, we presented our preliminary results on the study of the charm quark mass 
value within the framework of QCD Sum Rules at 0{al). Our hnal numbers are still 
preliminary and are not reported in the present discussion. We proht from the discussion 
on our results with our experimental colleagues present in this workshop. We are thankful 
to the organizers for providing such an excellent environment for that kind of interactions. 
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2.7 Current status of carlomat_3.0, an automatic tool for low en¬ 
ergetic electron-positron annihilation into hadrons 

K. Kolodziej 

Institute of Physics, University of Silesia, ul. Uniwersytecka 4, PL-40 007 Katowice, Poland 


In the energy range below the J/tjj threshold, processes of e’''e“-annihilation to hadrons 
cannot be described in the framework of perturbative QCD. The scalar electrodynamics 
(sQED) which was implemented in carlomat_2.0 [T] does not describe them in a satisfactory 
manner either. The theoretical frameworks usually used for the description of such processes 
are the Resonance Chiral Theory (RyT) or Hidden Local Symmetry (HLS) model which 
were proven to be equivalent in this context [2]. Both the RyT and HLS model involve, 
among others, the photon-vector meson mixing and a substantial number of vertices of a 
rather complicated Lorentz tensor structure that is present neither in the Standard Model 
(SM) nor sQED. Even at low energies, the final state of e’''e“ —)■ hadrons may consist of 
several particles, such as pions, kaons, or nucleons with one or more photons, or light fermion 
pairs such as or receive contribution of a big number of the Feynman diagrams. 

A new version of a program carlomat [3], tagged with index 3.0 |1], is just dedicated to the 
description of such multiparticle processes involving hadrons in a fully automatic way. 

Implementation of the photon-vector meson mixing required substantial modification 
of the code generation part of carlomat. Further changes in that part of the program 
were needed in order to incorporate calls to new subroutines for computation of the helicity 
amplitudes of the building blocks and complete Feynman diagrams which contain new inter¬ 
action vertices and mixing terms. Moreover, a number of new subroutines for computation 
of the helicity amplitudes involving the Feynman interaction vertices of the RyT or HLS 
model [3] with the Lorentz tensor structures, which are not present either in the SM or in 
the effective models implemented in the former version of the program, were written and 
many subroutines have been modified in order to incorporate the g^-dependent couplings 
and vector meson widths. 

In order to give the user a better control of the implemented effective models, a number 
of new options have been introduced in the Monte Carlo computation part of the program. 
They include a possibility of choosing different formulae for the g^-dependent couplings, 
where the four momentum transfer g is determined automatically at the stage of code 
generation from the four momentum conservation in the corresponding interaction vertex, 
or for the s-dependent widths. It is also possible to modify weights and phases of each of the 
photon-vector meson mixing terms implemented in the program, which should help to find 
out the dominant production mechanisms of different hadronic channels. Another important 
option offers a possibility to test the electromagnetic gauge invariance for processes involving 
one ore more photons. 

carlomat_3.0 can be downloaded from CPC Program Library or from the web page: 
http://kk.us.edu.pi/carlomat.html. 

Acknowledgement: This project was supported in part with hnancial resources of the 
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The author is grateful to Fred Jegerlehner for providing the Feynman rules of the HLS model. 
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2.8 A first glance towards a dispersive formalism for 7*7 —> tttt 


P. Masjuan 

PRISMA Cluster of Excellence, Institut fiir Kernpliysik and Helmholtz Institut Mainz, 
Johannes Gutenberg-Universitat Mainz, D-55099 Mainz, Germany 


We present a hrst step towards a dispersive formalism for the 7*7 process. This 

process is not only interesting by its own, as it will be measured for both one- and two-virtual 
photons in BESIII and Belle with high precision and encodes at once several interesting 
theoretical aspects (gauge invariance, hnal-state interactions, form factors), but also for 
its potential relation to the hadronic light-by-light scattering (see and references 

therein) and the generalized pion polarizabilities. 

Our goal is to provide a friendly useful MonteCarlo parameterization based on dispersion 
relations gl El El El E] while keeping its essential ingredients [9]. This demands identifying 
the crucial pieces of information that allow for a reliable, albeit not complete, description of 
the current mm and forthcoming data at BESIII. The analysis of the recent works iiEii, 
and taking into account only the tttt channel (neglecting then any inelasticity up to almost 
1 GeV), suggests to neglect the KK channel, while keeping left-hand cut contributions 
beyond the one-pion exchange in a manageable way na. 

With this ingredients at hand, we proceed to describe an unsubtracted dispersion re¬ 
lation (DR) as in Ref. [5] but with both S- and D-waves unitarized using, thanks to the 
Fermi-Watson theorem, the corresponding tttt phase-shifts, solutions of which are taken from 
Ref. [13] for both isospin J = 0,2 (the only accessible with two incoming photons). Since, 
however, we want to consider the elastic tttt channel only, we should modify the phase shift 
around the KK threshold for emulating, without employing a coupled channel formalism, 
the amplitude’s phase. The coupled channel is known lai but difficult to implement with 
virtual photons [12] , and hard to simplify for a user friendly MonteGarlo generator. For that 
we use the proposal of Ref. [12] to dehne a piecewise function. The high energy behavior 
of the phase-shift, beyond the parameterization provided in [T3], is matched to the Regge 
behavior [ 8 ] using the latest Regge analysis of the PDG data HI. 

Virtual photons imply extending the DR for 77 —)■ tttt, which has two independent 
helicity amplitudes, to three (one-virtual photon) and hve (two-virtual photons) independent 
amplitudes. For the one-virtual photon discussed in this talk, the new amplitude yields an 
enhancement on the cross section around the threshold of the same order as the transversal 
amplitudes. Such effect helps to compensate the decrease of the cross section due to the flux 
factor suppression in presence of photon virtualities together with the 7 — form factor 
(and the soft-photon limits [T 2 ]b 

As we just said, as an ingredient we need the coupling of the virtual photon with the 
pion (given by the vector form factor). We use a data driven parameterization specially 
suited for low energies from Ref. Ha . Around the /2(1275) tensor resonance, which shows 
up in the the 77 —)■ tttt process around l.SGeV, the amplitude can be worked out with 
the help of a Breit-Wigner model for the /2 [5], or using the appropriate tttt phase-shifts 
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together with the left-hand-cuts of the amplitudes (i.e, vector and axial-vector contributions 
in the t—channel exchange). Even though we are concerned with the low-energy sector of 
the process, we also include the effects of the /2 in our formalism since they spread all over 
the cross section. Within the former approach, the /2 coupling to the virtual photon can 
be parameterized in terms of vector meson dominance ideas [16], using the rj and rj' form 
factors im. For the later, the coupling of the virtual vector with the pion and the (left- 
hand) vector/axial-vector, in terms of meson dominance [TBl [T2]. For a photon virtuality 
around = 0.5 GeV^, the typical intensity of the integrated cross section is of around 30 
nb [9], measurable at BESIII. 
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2.9 Current status of two and three pion decay modes within 
TAUOLA 


O. Shekhovtsova 

Institute of Nuclear Physics PAN ul. Radzikowskiego 152 31-342 Krakow, Poland 
Kharkov Institute of Physics and Technology 61108, Akademicheskaya,!, Kharkov, Ukraine 


Since 90’ years Tauola is the main Monte Carlo generator to simulate tan-lepton decays 
[T]. It has been used by the collaborations Aleph [2], Cleo [5], at both B-factories (BaBar 
[1] and Belle [5]) as well at LHC [B] experiments for tan decay data analysis. In view of the 
forthcoming Belle-II project [7] it is important to revise the Tauola context in details. In 
this notes we discuss the status of the predominant hadronic tan-lepton decay modes: two 
pion [Br ~ 25.52%) and three pion [Br ~ 18.67%) decay channels. 

In the general case, a hadronic current of two-meson tan-lepton decay mode depends on 
both vector and scalar form factors. However, in the isospin symmetry limit, m^- = mj^o, 
the scalar form factor vanishes for two-pion decay mode and the current is described by 
the vector form-factor only. Currently Tauola has four parameterizations for the vector 
form-factor: 

• Kuhn-Santamaria (KS) parametrization [ 8 ]; 

• Gounaris-Sakurai (GS) parametrization used by Belle [5], Aleph and Cleo collabora¬ 
tions; 

• parametrization based on the Resonance Chiral Lagrangian (RChL) [9]; 

• combined parametrization (combRChL) that applies dispersive approximation at low 
energy and modified RChL result at high energy HD]. 

In all cases, except for the RChL parametrization, the pion form factor is given by interfering 
amplitudes from the known isovector meson resonances p(770), p'(1450) and p"(1700) with 
relative strengths 1, (3 and 7 . Although one could expect from the quark model that jS and 
7 are real, their phases are left free in the fits. In the case of the RChL parametrization only 
the p(700) and p'(1450) contributions are included, with the relative p' strength being a real 
parameter. For the energy-dependence of p(770) two-pion and two kaon loop contributions 
are included for RChL and combRChl parameterizations, whereas in the case of KS and 
GS the p width is approximated by the two pion loops only. The p' and p” widths include 
two-pion loops only for all parameterizations. 

Results of the fit to the Belle data [5] are presented in Fig. 1. The best fit is within 
the GS pion form-factor and the worst one is within the RChL one. To check the influence 
of the p"(1700) on the RChL result, the p” resonance has been included in the same way 
as it was done for p'; however, this inclusion has not improved the result. Therefore, we 
conclude that missing loop contributions could be responsible for the disagreement and that 
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Figure 1; The pion form factor fit to Belle data [5]: the GS parametrization (left panel), 
the RChL parametrization (at centre), right panel is with the combRChL parametrization 
(right panel). 


the phases of [5 and 7 might mimic the missing loop contribution. This idea will be checked 
by adding, hrst, a four-pion loop contributions. 

When the r lepton decays into three hadrons and a neutrino, the predominant decay 
mode involves three pions. In the general case Lorentz invariance determines the decompo¬ 
sition of the hadronic current for a three hadron hnal state in terms of hve Lorentz invariant 
structures [H] multiplied by hadronic form factors. Among the three hadronic form factors 
which correspond to the axial-vector part of the hadronic tensor, only two are independent. 
In Tauola the following three pion form factors are available 

• CPC version [ 1 ], that includes only the dominant ai pvr mechanism production 

• Cleo parametrization with equal currents for and 7 r“ 7 r“ 7 r’'' modes based on [3] 

• Cleo parametrization with the current [3] and the current from the 

unpublished Cleo analysis [12] (the currents coincide for vector and axial-vector inter¬ 
mediate states and are not equal when the scalar and tensor resonances are included) 

• modihed RChL parametrization |T3| 

The Cleo parametrization results have not yet been compared with the BaBar prelim¬ 
inary data [1] and it will be a task of future work. Whereas two- and three-pion spectra 
calculated on the base of the modihed RChL parametrization have been htted to the BaBar 
preliminary results. Discrepancy between theoretical spectra and experimental data can be 
explained by missing resonances in the model, such as the axial-vector resonance a'^(1600), 
the scalar resonance /o(980) and the tensor resonance /2(1270). Inclusion of these resonances 
in the RChL framework will be a future task. 

Comparison of the 7 r“ 7 r“ 7 r’'' current in the framework of the modihed RChL with the 
ChPT result has demonstrated that the scalar resonance contribution has to be corrected 
to reproduce the low energy ChPT limit. The corresponding calculation is under work. 
The same type of comparison has to be done for both version of the Cleo currents 

described above. 
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2.10 A new approach to evaluate the leading hadronic corrections 
to the muon g-2 
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The long-standing discrepancy between experiment and the Standard Model (SM) pre¬ 
diction of a^, the muon anomalous magnetic moment, has kept the hadronic corrections 
under close scrutiny for several years B121 El S]. In fact, the hadronic uncertainty domi¬ 
nates that of the SM value and is comparable with the experimental one. When the new 
results from the g-2 experiments at Fermilab and J-PARC will reach the unprecedented 
precision of 0.14 parts per million (or better) B El Eji fhe uncertainty of the hadronic 
corrections will become the main limitation of this formidable test of the SM. 

An intense research program is under way to improve the evaluation of the leading 
order (LO) hadronic contribution to a^, due to the hadronic vacuum polarization correction 
to the one-loop diagram as well as the next-to-leading order (NLO) hadronic one. 

The latter is further divided into the 0{a^) contribution of diagrams containing hadronic 
vacuum polarization insertions [10], and the leading hadronic light-by-light term, also of 
0{a^) BlIIllIl]. Very recently, even the next-to-next-to leading order (NNLO) hadronic 
contributions have been studied: insertions of hadronic vacuum polarizations were computed 
in [13], while hadronic light-by-light corrections have been estimated in [14] . 

The evaluation of the hadronic LO contribution a™° involves long-distance QCD for 
which perturbation theory cannot be employed. However, using analyticity and unitarity, it 
was shown long ago that this term can be computed via a dispersion integral using the cross 
section for low-energy hadronic e+e“ annihilation [T3]. At low energy this cross-section is 
highly fluctuating due to resonances and particle production threshold effects. 

An alternative determination of can be obtained measuring the effective electro¬ 
magnetic coupling in the space-like region extracted from Bhabha e~^e~) scattering 

data [16]. As vacuum polarization in the space-like region is a smooth function of the squared 
momentum transfer, the accuracy of its determination is only limited by the statistics and 
by the control of the systematics of the experiment. Also, as at flavor factories the Bhabha 
cross section is strongly enhanced in the forward region, we will argue that a space-like de¬ 
termination of may not be limited by statistics and, although challenging, may become 
competitive with standard results obtained with the dispersive approach via time-like data. 
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